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Abstract. We investigate the sample paths regularity of operator scaling a-stable random fields. 
Such fields were introduced in [6] as anisotropic generalizations of self-similar fields and satisfy the 

". scaling property {X{c^x);x G R"^} ^^=^ {c"X{x)-x £ W^} where E is a, dxd real matrix and H > 0. 

' In the case of harmonizable operator scaling random fields, the sample paths are locally Holderian 

I and their Holder regularity is characterized by the eigen decomposition of R"^ with respect to E. In 

■ particular, the directional Holder regularity may vary and is given by the eigenvalues of E. In the 

^ ' case of moving average operator scaling random a-stable random fields, with a £ (0, 2) and d > 2, 

^ , the sample paths are almost surely discontinous. 

CN ' 1- Introduction 



In this paper we consider operator scaling stable random fields as introduced in [6j. More precisely, 
for E a real dxd matrix with positive real parts of the eigenvalues, a scalar valued random field 
is called operator scaling for E and H > ii for every c > 



(1) {X{c^xy,x eR'^}^=^'' {c^X{x);x eR'^}, 

where ^'^= ^ means equality of finite dimensional distributions and as usual = exp(Slogc) with 
exp(A) = X^^Q 4r is the matrix exponential. Let us remark that up to consider the matrix E/H, 
^ ' we may assume, without loss of generality, that H = 1. 

. These fields can be seen as anisotropic generalizations of self-similar random fields. Let us recall 

^ I that a scalar valued random field {X{x))^^^d is said to be H -self-similar with G M if 



O : {X{cx)- X £ W^} ^ ='^^ {c"X{xy, x G M^} 

, for every c > 0. Then a iJ-self-similar field is also an operator scaling field for the matrix E = Id/H, 
I where Id is the identity matrix of size dxd. Self-similar random fields are used in various applications 
■ ■ such as internet traffic modelling [29], ground water modelling and mathematical finance, just to 

mention a few. Various examples can be found for instance in the books [THKHEZ]- A very important 
class of such fields are given by Gaussian random fields and especially by the fractional Brownian 
field Bh, where H € (0, 1) is the so-called Hurst parameter. The random field Bh is -ff-self-similar 

^ ■ and has stationary increments, e.g. {Bh{x + /i) — BH{h);x G M'^} ^'^=^ {Bh{x);x G M"'} for any 
/i G M'^. It is an isotropic generalization of the famous fractional Brownian motion, implicitely 
introduced in [12] and defined in [19] . 

Nevertheless, Gaussian random fields are not convenient for some heavy tails phenomena mod- 
elling. For this purpose, a-stable random fields have been introduced. A scalar valued random field 
{X{x);x G M"^} is said to be symmetric a-stable (SaS), for a G (0,2), if any linear combination 

n 

J2 O'kXixf^) is SaS. We address to the book [27] for a well understanding of such fields. The self- 

k=l 

similar a-stable fields with stationary increments have been extensively used to propose alternative 
to Gaussian modelling (see [231 122] for instance) and are isotropic. 

However, certain applications (see, e.g., [3 [7] and references therein) require that the random 
field is anisotropic and satisfies a scaling relation. This scaling relation should have different Hurst 
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indices in different directions and these directions should not necessarily be orthogonal. To reach 
this goal the authors of [6] propose a new class of random fields with an anisotropic behavior driven 
hy a d X d matrix E. More precisely, they introduce a-stable random fields which have stationary 
increments and satisfy the operator scaling property Two different classes of such fields are 
defined and analyzed, using a moving average representation as well as an harmonizable one. In 
the Gaussian case a = 2, according to [6] there exist modifications of these fields which are almost 
surely Holder-continuous of certain indices. We give similar results here in the stable case a G (0, 2) 
for harmonizable operator scaling stable random fields. Actually, we obtain their critical global and 
directional Holder exponents, which are given by the eigenvalues of E. In general, such fields are 
anisotropic and their sample paths properties varies with the direction. In particular, in the case 
where E is diagonalizable, for any eigenvector 9j associated with the real eigenvalue \j, harmonizable 
operator scaling stable random field admits Hj = 1/Xj as critical Holder exponent in direction 9j. 
Let us point out that we establish an accurrate upper bound for the modulus of continuity. Such 
upper bound has already been given in the case of real harmonizable fractional stable motions, which 
are self-similar, in [T3] and in the case of some Gaussian random processes in |13] . Then, in this 
paper, we generalize these results to any dimension d and any harmonizable operator scaling stable 
fields. From this upper bound, we can deduce the Holder sample paths regularity. Let us point out 
that we also obtain this upper bound in the case of Gaussian operator scaling random field, which 
improves the sample paths properties establishes in [6]. 

Furthermore, whereas in the Gaussian case a = 2, moving average and harmonizable fields have 
the same kind of regularity properties, this is no more true in the case a G (0,2). In particular, 
we show that for d > 2, a moving average operator scaling stable random field does not admit 
any continuous modification. Remark that if d = 1, the sample paths regularity properties are 
already known since the processes studied are self-similar moving average stable processes, see for 
example (271 H [28] • 

One of the main tool for the study of sample paths of operator scaling random fields is the change 
of polar coordinates with respect to the matrix E introduced in [22] and used in [6]. For X a Gaussian 
operator scaling random field with stationary increments, using ([T]), we can write the variogramme 
of X as 

v\x) = E{X^{x)) = TE{xf''E{X\£E{x))), 

where te{x) is the radial part of x with respect to E and £e{x) is its polar part. Therefore, in 
the Gaussian case, the sample paths regularity depends on the behavior of the polar coordinates 
{te{x),Ie{x)) around x = 0. Such property also holds in the stable case a S (0,2). The Holder 
regularity of the sample paths follows from estimates of te{x) compared to ||2;||. These estimates are 
given in Section [3] and their proof are postponed to the Appendix. 

Furthermore, to study the sample paths in the stable case, the other main tool we use is a series 
representation of harmonizable operator scaling a-stable random fields. Such representations in 
series of infinitely divisible laws have been studied in [T71 [TH [23 [M] • As in [H] , our study is based 
on a LePage series representation. Actually, the main idea is to choose a representation which is a 
conditionnally Gaussian series. 

In Section [2l we recall the definition of harmonizable operator scaling random fields. Then, Sec- 
tions [3] and [H are devoted to the main tools we need for the study of the sample paths of these 
fields. More precisely. Section [3] deals with the polar coordinates with respect to a matrix E and 
Section [H gives the LePage series representation. In Section [5l the sample paths properties and the 
Hausdorff dimension of the graph are given. Section [6] is concerned with moving average operator 
scaling random fields. 



2. Harmonizable representation 
Let us recall the definition of harmonizable operator scaling stable random fields, given by [0]. 
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Let E he a real d x d matrix. Let Ai, . . . , be the complex eigenvalues of E and aj = 5?(Aj) for 
each j = 1, . . . ,d. We assume that 

(2) min a,- > 1. 

i<j<rf 

Let ip : M!^ — > [0, oo) be a continuous, i^^-homogeneous function, which means according to Defini- 
tion 2.6 of [6] that 

ip{c^ x) = cij){x) for all c > 0. 

We assume moreover that ^'(2;) 7^ for 2; 7^ 0. Such functions were studied in detail in [22], Chapter 5 
and various examples are given in Theorem 2.11 and Corollary 2.12 of [B]. 

Let < a < 2 and Woi{d^) be a complex isotropic a-stable random measure on with Lebesgue 
control measure (see [27j p. 281). If a = 2, Wa{dS,) is a complex isotropic Gaussian random measure. 
Let q = trace {E). 

Definition 2.1. Since ([2]) is fulfilled, the random field 

(3) X^{x) = K / (e*<^'«> - l)^(e)-i-'^/" W^{di) ,x G 

is well defined and called harmonizable operator scaling stable random field. 

From Theorem 4.1 and Corollary 4.2 of [6], is stochastically continuous, has stationary incre- 
ments and satisfies the following operator scaling property 

(4) Vc > 0, {X«(c^x);x G R'^^ ^ ='^^ {cX„(x);x G R'^}. 
For notational sake of simplicity we denote the kernel function by 

(5) /(;^,^) = fe*(-,0_i)^(^)-W". 



Let us recall that f{x,-) G L"(M'^) for any x G M.'^, which is a necessary and sufficient condition 
for Xa to be defined. 

Now, let us give some examples of such random fields. 

Example 2.1. Let Id be the identity matrix of size d x d, E = Id/ H (with Q < H <1) and ip[x) = 
1 1 2; 1 1 with II • II the Euclidean norm. Then the random field defined by ([3]) is a real harmonizable stable 
random field (see [27j for details on such fields). In this case, X^ is self-similar with exponent H, i.e. 

Vc > 0, {X„(cx);x G R'^} ^=^^ {c^X„(x);x G M"^}. 

Let us quote that, if a = 2, X^ is a fractional Brownian field and its critical Holder exponent is given 
by its Hurst index H (see Theorem 8.3.2 of [2] for instance). 

Example 2.2. Assume that E is diagonalizable. Then, all the eigenvalues of E are real given 
by (aj)i<j<^- Let {dj)i<j<d be a basis of some corresponding eigenvectors and consider the function 
defined by 

The function is clearly continuous and non negative on R'^. Moreover, since (c^*x, 6j) = (x, c^9j) = 
c"^ (x, 9j), it is also i^^-homogeneous. Finally, since {Oj)i<j<d is a basis of R'^ we have that tp^x) = 
if and only if x = 0. Then we can define X^ by ^ and in this case the operator scaling property ^ 
implies that 

Vj = 1, . . . , d, Vc> 0, {Xa{ct9j);t G R} ^■^= ^ |c^/''^X«(t6'j); t G m}. 

The random field X^ is self-similar with Hurst index Hj = l/aj in the direction 9j. In particular, in 
the Gaussian case (a = 2), {X2{t9j))^^^ is a fractional Brownian motion with Hurst index Hj. Then, 
in this case, its critical Holder exponent is equal to Hj. 
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The main tool in the study of operator scaling random fields is the change of coordinates in a kind 
of polar coordinates with respect to the matrix E. Then, before we study the sample paths of X^., 
we recall in the next section the definition of these coordinates and give some estimates of the radial 
part. 

3. Polar coordinates 

According to Chapter 6 of [22]) since E is a real d x d matrix with positive real parts of the 
eigenvalues there exists a norm || • H^; on M'^ such that for the unit sphere Se = {x E M'^ : \\x\\e = 1} 
the map 

^E : (0, oo) xSe M*^ \ {0} 

{r,e) I — > r^e 

is a homeomorphism. Hence we can write any x S M'^\{0} uniquely as x = te{x)^£e{x) for some 
radial part te{x) > and some direction iE{x) S Se such that x i— > te{x) and x i-^ (.e{x) are 
continuous. Observe that S^; = {x G IR'^ : te{x) = 1} is compact and set 

(6) iriE = min||x|| and Me = max||x||. 

Se Se 

We know that te{x) ^ oo as x — > oo and te{x) — > as x — > 0. Hence we can extend te{-) 
continuously by setting te{0) = 0. 

Let us recall the formula of integration in polar coordinates established in [6]. 

Proposition 3.1. There exists a unique finite Radon measure ue on Se such that for all f G 
L^(E.'^, dx) we have 



/ f{x)dx= / / f{r'^e)aE{d9)r^~^dr. 
jR'i Jo J Se 



In the Gaussian case (a = 2), the variogramme of X2 can be rewritten as follows 



TE{xf^[X2{^E{x)y 



v'^{x) 

using the operator scaling property Then, the Holder regularity of X2 follows from estimates 
of te{x) compared to ||x|| around x = 0, e.g. the Holder regularity of te around 0, see [6]. Then, 
in order to get an upper bound for the modulus of continuity (for any a), we need precise estimates 
of te{x). 

As done in [21] for the study of operator-self-similar Gaussian random fields we use the Jordan 
decomposition of the matrix E to get estimates of te- From the Jordan decomposition's theorem 
(see [To] p. 129 for instance), there exists a real invertible d x d matrix P such that D = P~^EP 
is of the real canonical form, which means that D is composed of diagonal blocks which are either 
Jordan cell matrix of the form 

(\ 1 ... 0^ 

A 



V ... 

with A a real eigenvalue of E or blocks of the form 



1 

A / 







(7) 



h 
A 




h 



\ 





h 

A / 



with A 



and h 



\ 

where the complex numbers azizib (6 7^ 0) are complex conjugated eigenvalues of E. 
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Let us denote by || • || the subordinated norm of the Euchdean norm on the matrix space. Precise 
estimates of te follow from the next lemma. 

Lemma 3.2. Let J he either a Jordan cell matrix of size I or a block of the form ([7]) of size 21 
associated with the eigenvalue X = a + ib. Then, for any t G (0, e~^] U [e, +oo) 

t" < lli^ll < V2/et"|loet|'-^ 



Proof, see the Appendix. 



□ 



Let us be more precise on the Jordan decomposition of E. Let us recall that the eigenvalues of 
E are denoted hy Xj, j = 1 . . . d and that aj = 5?(Aj) > 1 for j = 1, . . . , d. There exist Ji, . . . , Jm, 
where each Jj is either a Jordan cell matrix or a block of the form ([7]) , and P a real d x d invertible 
matrix such that 

/ Ji ... \ 

J2 : 



E = P 



p-\ 



: ••. ■•. 

V ... Jm / 

We can assume that each Jj is associated with the eigenvalue Xj of E and that 

1 <«!<••• < am- 



We also set Hj = a^ ^ and have 



(8) < i7„ < ••• < -H^i < 1. 

If Xj G M, Jj is a Jordan cell matrix of size Ij = Ij G N\{0}. If Xj G 
of size Ij = 21 j G 2N\{0}. Then for any t > 0, 

/ t-^i ... \ 

: 



, Jj is a block of the form ([7|) 



t^ = P 



V 



p- 



t-^™ / 



We denote by (ei,...,erf) the canonical basis of M'^ and set fj = Psj for every j = 1,, 
Hence, (/i, . . . , fa) is a basis of W^. For all j = 1, . . . , m, let 

(9) Wj = Vect (fk ■,Y.k + l<k< Y,l]\ . 



1=1 



1=1 



Then, each Wj is a £^-invariant set and E = ©T^i Wj. 



The following result gives bounds on the growth rate of te{x) in terms of the real parts of the 
eigenvalues of E. 

Proposition 3.3. For any r G (0, 1) there exist ci,C2 > such that for every 1 < < j < m, 

Ci||x||^^0 |log ||x|||"^^^O.i-l)^io < Te{x) < C2\\x\\^^ \\og\\x\\\^^'0'0-'^)^i 



holds for any x G with \\x\\ < r and pj^j = max 

■"^ ' jo<k<j 

Proof. See the Appendix. 

Then, we easily deduce the following corollary. 
Corollary 3.4. For any r G (0, 1) there exist ci,C2 > such that for any x G Wj with \\x\\ < r 

(10) cillxll^^- llogllxlir^'^"^)^^' < te{x) < C2\\x\\"' |log||2;|||('^"^)^^' 
and for any x G M"^ with \\x\\ < r 

(11) cillxll^i |log 11x11 r^'~^^^^ < te{x) < C2||x||^™ |log||x|||('"^)^'" 



□ 
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where I = max Ij. 

l<j<m 



Therefore we have precise estimates for the Holder regularity of the radial part. Let us remark that 
we improve the first statement of Lemma 2.1 of [6] and that the second one can also be improved in 
a similar way. From these estimates we deduce the Holder regularity of in Section [5l As already 
mentionned, the study of the sample paths is based on a series representation. Then, before we state 
regularity properties, it remains to give the LePage series representation of harmonizable operator 
scaling random fields Xa defined by ([3|). 



An overview on series representations of infinitely divisible random variable without Gaussian part 
can be found for example in [24^ [26] and references therein. In particular, LePage series represen- 
tation ([TTldB]) have been used in [I5l [Hj to study the sample paths regularity of some self-similar 
a-stable random motions with a £ (0,2). Here, this representation is also the main representation 
we use in the case a G (0, 2). Actually, in the Gaussian case q = 2, such representation does not hold. 

Let us now introduce some notations that will be used throughout the paper. Let /i be an arbitrary 
probability measure equivalent to the Lebesgue measure on M"^ and let m be its Radon-Nikodym 
derivative that is /i((i^) = m(^)d^. 

Notation Let (T'«)„>i, (5'n)„>i and iS,n)n>i be independent sequences of random variables. 

• Tn is the nth arrival time of a Poisson process with intensity 1. 

• {gn)n>i is a sequence of i.i.d. isotropic complex random variables so that gn == e^^gn for any 
6* G M. We also assume that < Ed^-nl") < +oo. 

is a sequence of i.i.d. random variables with common law n{d^) — m{^)d^. 

According to Chapter 3 and Chapter 4 of |27j . stochastic integrals with respect to an a-stable 
random measure A can be represented as a LePage series as soon as the control measure of A is a 
finite measure. The next proposition generalizes this to Wa whose control measure is the Lebesgue 
measure. It is a consequence of Lemma 4.1 of [15], which is a correction of Lemma 1.4 of [20]. This 
proposition can also be deduced from [241 125] . concerned with series representations of stochastic 
integrals with respect to infinitely divisible random measures. 

Proposition 4.1. Assume that a G (0,2). Then, for every complex valued function h G L"(M'^), the 



4. Representation as a LePage series 



series 



(12) 




n=l 



converges almost surely. Furthermore, 




with 




1/q 



Remark 4.1. According to Proposition \4-l\ taking a G (0,2), the random measure 



^a{dO = CaY^T-'/^miUr'^'^grAnidC) 



n=l 



is a complex isotropic a-stable random measure. 
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Proof. Let K = m{^n)~^/"h{C )gn- Then, n > 1, are i.i.d. isotropic complex random variables. 
By Lemma 4.1 in [1^, converges almost surely and 



Vz G C, E(^exp (i'Sl(zY^'^jj = exp (-cj"|z|"; 

with 

Jo ^ 

Since gi is invariant by rotation and independent with ^i, 

nmviT) = Efm(6)">(6)r)iE(i??(5i)r) = nmgiT) I mrac. 



Moreover, by definition of an isotropic a-stable random stable measure (see 

yz eC,E(^exp(^m(^z j h{i)M{dC)^y^ = exp(-c2(/i)|zr) 

with c"(/i) = ^— j \cos {x)\'^ dx^ J^^ \h{S,)\°'dS,- Therefore, we have 

where Ca is defined by (|13p . which concludes the proof. □ 

From the previous proposition, we deduce the following statement which is the main series repre- 
sentation we use in our investigation. 

Proposition 4.2. Let a G (0, 2). For every x G W^, the series 

(14) y„(x) = Ca^(^T-'/''m{Uy'^''f{x,^n)g}j , 

where Ca is defined by (jl3|) . converges almost surely. Furthermore, 

|y„(x) : X G R'^} ^=^^ {Xa{x) : x G M"^}. 

Proof. From Proposition 14.11 for any x G M'^, the convergence of the series follows from the fact that 
f{x, •) G L"(IR'^). The equality in distribution between Xa and Ya is obtained by linearity of both 
fields. □ 

Using LePage representation (jl4p of and the estimates given in Section [3l we give an upper 
bound for the modulus of continuity of Xa and obtain the critical Holder regularity of its sample 
paths in the next section. 

5. Holder regularity and Hausdorff dimension 

Throughout this section we fix a compact set of M*^ and investigate the Holder regularity of the 
sample paths of X^ on K, with X^ the harmonizable operator scaling stable random field defined 

by m- 

Let us recall that for the Gaussian case a = 2, according to Theorem 5.4 of [6j, the Holder regularity 
of X2 depends on the subspaces (^j)i<j<m defined by ^ and associated to the eigenvalues of E. 
More precisely. Theorem 5.4 of [6] implies that, when restricted to the subspace Wj, the Gaussian 
random field {X2{x);x G Wj} admits Hj as critical Holder exponent. This follows from the fact 
that the regularity of X2 is determined by the regularity of te around 0, which is given by Hj 
according to (|1U|) . Here, we give an upper bound for the modulus of continuity of X^ in the general 
case a G (0, 2]. Then we prove that the critical Holder exponents are the same than in the Gaussian 
case a = 2. Let us state our main result. 
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Theorem 5.1. Let a G (0,2). There exists a modification X* of on K such that 
(15) lim sup KM-X:fa)| ^ „ ^ ^ 

Te{i - !/)|logTE(l - >,)|"»+'/2+' 

for any e > 0. 

This resuh was proved in the case of harmonizable self-similar processes in [T3], e.g. in the case 
of Example 12.11 with d = 1. The main idea is to use the LePage series representation (jl4p where Qn, 
n > 1, are Gaussian complex isotropic random variables. Furthermore, it remains to choose the 
density distribution m of the In [H], the authors choose 



iei|log|C|| 

where c,j > 0. A straightforward generalization in higher dimension d leads to choose 

m(0 = — — 7^, ^ G M'^\{0}. 

Remark that in this case (e.g. Example 12. Ih the matrix E = Id/H = and that we can choose 
II ■ He* = II ■ II- Then, using classical polar coordinates, we obtain that for x / 0, 

TEtix) = \\x\\^ and Ie^-{x) = 71 — n" 

||x|| 

and therefore 



rs*(?niogrE*(Ol'+'' 

since q = trace(£') = d/H. Then we choose this density in the general case. The advantage is that m 
only depends on the radial part r^;*. 

Proof of Theorem \5.1[ We can assume without loss of generality that K = [0,1]'^. According to 
Proposition 14.21 for every x G R'^ 

y„(x) = C„5R|^^^r-^/"?n(e„)-'/"/(x,en)5n^ 

converges almost surely and Ya '''= Xa- As already mentionned, we assume that gn, n > 1 are 
Gaussian complex isotropic random variables and that the density distribution of is defined by 

TEtiO \^OgTEt{0\ 

where 77 > and c,, is chosen such that J^^ m{S,)d^ = 1. In particular, conditionally to (T„, Ya{x) 
is a real-valued Gaussian random variable. 

As in the proof of Kolmogorov-Centsov Theorem (see |llj). we first prove that almost surely 
for TE{xk — Xk') small enough, 

\Ya{Xk) - Ya{xk')\ < CTE{Xk - Xk>)\log TE{xk - Xfc/) | . 

where (x^)^ is some countable dense sequence of K. Then, Xq. satisfies the same property. 
Finally, we give a modification X* of X^ for which (jlSp holds. In the first step, we construct the 
sequence (x^)^ we use. 

Step 1. For k £ N\{0} let us choose i^k £ I^\{0} the smallest integer such that 

C2d^'"/^2-'^''-^™|z^fclog2|('~^^^'" < 2"*^, 

where C2 and / are given by Corollary 13.41 for r = 1/2. Remark that by definition, i/k < J^fc+i- Up to 
change C2 in Proposition 13.41 we can assume that 
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which impHes that Vk> \ for every k. Then, since 1 < i^fc — 1 < z^fc, 
and we have 

(16) 2~^{2^~^^ c^^ < (^-""^ {uk\og2)^"^Y^ < 2~^[Vdy""^ ci^^ . 

Let us remark that {i^k)k>i increasing sequence and then that Vk>k for every k. Furthermore, 
taking the logarithm of (I16p . one easily proves that 

k 

hm — = Hm-, 

which imphes = 0{k). 

For every k G N\{0} and j = (ji, . . . , j^) G Z'' we set 

Xkj = ^,Vk = [xk,j ■.jeZ'^n [0, 2^^]'^} and iVfe = cardPk = (2"^ + 1)'^. 

Then 2?^ is a 2~'^ net of K for r^; in the sense that for any x K one can find x^j S such that 
TE{x—Xkj) < 2"^^. Actually, by Proposition l3.4l it is sufficient to choose j such that ji < 2'^'=Xi < 
for 1 < i < d. 

+ 00 

Let us remark that the sequence (T^k)k increasing and set = |^ T^k- 

k=l 

Step 2. Almost surely, for any x,y £T> 

where Ca is defined by (fT3]) . Since the sequences (Cn.)n and (s'n)n are independent and 

is a sequence of i.i.d. Gaussian complex isotropic random variables 

+ 00 

R{x,y) = 5;X-i/"m(e„)-^/"(/(x,e„) - f{y,Cn))9n 

n=l 

is a Gaussian isotropic complex random variable conditionally to {Tn,(,n)n- Remark that 1^0(3;) — 
Ya{y) = CcM{R{x-,y)) almost surely. Therefore, conditionally to {Tn,(,n)ni Yai^x) — Ya{y) is a real 
centered Gaussian random variable with variance 

v\{x,y)\(T^,Cn)J = ^E{\R{x,y)\'\{Tr.,Un) 

n=l 

since - /(y,Cn)| = \ f{x - y,in)\- 

We consider the set 

E^^j = {UJ : \Ya{Xk,i) - Ya{xkj)\ > v{{xk,i,Xk,j) I (Tn,(,n)n) ^irE{xk,i " Xkj))} , 

where, as in |13j . 

' 1 



(fit) = y2C^diog J, t > 

and > will be chosen later. Then, for every {k,i,j), 

We give in the following an upper bound of this probability for well chosen {k,i,j). Note that if Z 
is a real centered Gaussian random variable with variance 1, we have 

E^l^fe^ I (r„,^„)„^ = F{\Z\ > ^{TE{xk,i - Xk,j))) almost surely. 
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Let US choose 6 £ (0, 1) and set for k G N\{0}, 6k = 2-(i-'^)'= and 

Ik = {{i,3) G n [0,2*^^]^ : TE{xk,^ - Xk,j) < 5k]. 
For every G Ik, since 99 is a decreasing function 

F{\z\ > ipiTE{xk,^ - xk,j})) < n\z\ > ^m). 

Then, we recall that 

e-"'/2 

Vn > 0,F(Z > n) < . 

Therefore, for every k G N\{0} and G Ik, 



Hence, 



vr (/?(<5fc) v/2C7^d(l-5)A:log2 



y y Pfe,) <^=^^^=y2-(i-^)^^-'^cardIk. 

Since Vk = 0{k), the lower bounds of ()16p and Corollary 13.41 leads to 

cardjj eZ'^n [0,2'"= f : G 4} = O (^5^/^12'^'=/^™ yt^'^C"!)) , 
for any i G Z'^ H [0, 2'^'=]'^. Then one can find a finite constant C > such that 

which is finite for > and 6 small enough. By the Borel-Cantelli Lemma, almost surely 

there exists an integer k*{ijj) such that for every k > k*{u;), 

\Ya{x) - Ya{y)\ <v{{x,y)\ (Tn,^n)n) 'fi'^Eix - y)) 

as soon as x, y G Vk with TEix — y) < 6k- 



Step 3. As in [M] let us give an upper bound of 

C2 



\{x,y)\iTn,M = -f^{\9i\') T,T-'^''MCnr'/''\f{x-y,^n)\' 

n=l 

with respect to te{x — y). By definition of / 

v\{x,y) I {T^,in)n) < ^^{\gi\')a^ (TEix - y)) , 

where for all /i > 

+00 2 
a\h) = j;r-2/"m(^„)-2/"min [aIe \h'^\n\,2) m~'~''^", 

n=l 

with Me given by ([6]). For the sake of clearness we postpone the proof of the control of in 
Appendix and state it in the following lemma. 

Lemma 5.2. For any 7 G (0, 1) there exists a finite constant c > such that 

+00 

E(a2(/i)| (r„)J < cJ^r-2/"/i2|log/i|(^+'')(2/"-i) almost surely 

n=l 

as soon as h < 1 — 'j. 
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Following [13] let us denote 
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b{h) = /i|log/i|(^+'')/". 



Then by Lemma 15? 



HE 



(T„)„ < +00. 



Therefore by independence of (Tn)^ and (Cn)^, almost surely 



lim , „ , , , 



Up to change the Euclidean norm || • || by the equivalent norm || • \\^t defined in Lemma 6.1.5 of [22j 



the map h 



is increasing and so is a^. Also, one can conclude, as in [Tl], that almost surely 

lim -TTTTT- = 0. 



h-.o b'^{h) 

Therefore, up to change k* one can assume that for every k > k*{uj), for every x,y £ 2?^, 
(17) \Y^{x) - Y^iy)\ < y^TEix - y)\logTE{x - y)|(i+^)/"+i/2, 

as soon as te{x — y) < 5k- Let 



U n n <y2dC;Ts(x-y)|logrs(x-y)|(i+'')/"+i/2| 

n=l k>n Xjj/G'Dfe 

TE(^-V)<^k 

Since Xa and have the same finite dimensional margins P(r2*) = 1. 



Step 4. Let uj £ Vl* . By Step 3 there exists /c*(a;) > 1 such that satisfies ()17p for 
k > k*{oj), x,y GVk and Tsix - y) < 5k. 

Let us recall that by Lemma 2.2 of [6], there exists Ke > 1 such that for all x, y G 

te{x + y) <Ke {te{x) + TE{y)) ■ 

Let us denote F{h) = y2da^/i|log ^1^^+'')/"+^/^ and choose ko £ N such that 2''o5ko+i > 3if| 
and F is increasing on (O,^^,,]. Up to change k*{u>), we can assume that k*{uj) > k^. 



Let x,y £ V with x ^ y such that 3K^te{x — y) < 5k*(u))- Then there exists a unique k > k*{ui) 
such that 5k+i < 3K'^te{x — y) < 5k- Since x,y £ V, there exists n > k + I such that x,y £ Vn- 
Moreover, by Step 1, for j = k, . . . ,n — 1, we can choose x^^\y^^^ £ Dj such that 



te(x- x^^^i ) < 2~^ and te [y - y^^' ) <2~K 



By construction te [x^^^ — y^^^) < K"^ {te{x — y) + 22 . Let us point out that since k > /cq, 
2^5fc+i > 2''"5ko+i > 'iK'^- Therefore, one easily sees that 

TE (x^'''^ - y^'^A <^KlTE{x-y). 



Since 3K'^te{x — y) < 5k we obtain by Step 3 on the one hand that 
On the other hand we can write 



,ik) 



n-l 



j=k 
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with TE < 3K|;2~(-'+^) < 5j+i since j > k^. Moreover, note that x^^'^ G Vj C Pj+i 

and then by Step 3 

n-l 

+00 

where C = (j + §j < +cxd. Using similar computations for Xa{y) — Xa[y^^^), we 

j=0 

obtain that 

< {l + 2C)F{3KlTE{x-y)). 
Then one can find a constant C > such that for 3K'^te{x — y) < 
(18) \X^{x) - XUy)\ < Cte{x - y)\logTE{x - y)\i^+n)/-+y\ 

We now give a modification of X^. For 2; G D, we set 

X:{x){lo) = Xo^{x){u;). 

For X G if let x^'") G V such that lim„^+oo x^") = x. In view of (fTSj) . (X* (x^'")) (u;))^ is a Cauchy 
sequence and then converges. We set 



X:{x)iu;)= hrn X*(xW 

n— >+oo \ 



Remark that this limit does not depend on the choice of (x^"^). Moreover, since X^ is stochastically 
continuous, X* is a modification of X^. 
By continuity of te, we easily see that 

\X*{x){u;) - X*(x')(a;)| < Cte{x - y)|logTs(x - y)\i^+v)/a+i/2 

as soon as 3K'^te{x — y) < 5k* [ui)^ which concludes the proof. □ 

Following the same lines as the proof of Theorem 15.11 we obtain a similar result in the Gaussian 
case (a = 2) for more general fields. Let us remark that Ya is not defined for a = 2. However, in 
Step 2 of the proof, let us replace Ya hy X a. centered Gaussian random field and u^((x, y) \ {Tn, ^n)n) 
by the variance of X{x) — X{y) 



vH{x,y))=E[{X{x)-X{y)f). 



Furthermore let us replace Step 3 by the assumption that for some /3 G M and 6 > there exists a 
finite constant C > such that for x,y £ K with te{x — y) < 6 

(19) K{iX{x) - X(y))2) < CTEix - yfllogTEix - y)|^ 
Then Step 1, Step 2 and Step 4 yields the following proposition. 

Proposition 5.3. Let X be a centered Gaussian random field satisfying ([H]). There exists a modi- 
fication X* of X on K such that 

(20) hm sup ^ a.s. 

x,yeK TE{x-y)\logTE{x-y)\^/^+^^' 

\\x-y\\<S 

for any e > 0. 
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Let US point out that if X2 is an operator scaling Gaussian random field as defined in [6], then 

e((X2(x) - X2{y)f) = te{x - yf¥.(x2{lE{x - y)f 

and X2 satisfies (jl9p with /3 = by (5.2) of [6]. Therefore this result is more precise than one could 
expect from the Theorem 15. H replacing a by 2. 

Let us also mention that Marianne Clausel gives a different proof of a similar result for some 
Gaussian operator scaling random fields with stationary increments in [8]. 

Now, as in [6], we are looking for global and directional Holder critical exponents of the harmo- 
nizable stable random field Xa- These exponents have been introduced in [7] in the Gaussian realm 
but can be defined for any random field. Following Definition 5.1 of p], G (0, 1) is said to be the 
Holder critical exponent of a random field {X[x))^^-^d if there exists a modification X* of X such 
that for any s S {0,H), the sample paths of X* satisfy almost surely a uniform Holder condition of 
order s on i^, that is there exists a finite positive random variable A such that almost surely 

(21) \X* {x) - X* {y)\ < A\\x - yW" for all x, y G K 

while for any s G {H, 1), almost surely (|21|) fails. Note that the Holder critical exponent, if exists, 
is well defined since any continuous modification of X and X* are indistinguishable. Moreover, 
according to Definition 5.3 of [B] we say that X admits H{u) as directional regularity in direction 
u S S'^~^, with S"^~^ the Euclidean unit sphere, if the process {X{tu))^^^ admits H{u) as Holder 
critical exponent on K f] Mn. 

For all J = 1, . . . , 771 we set Kj = K f] ©'^^j^ Wk- Let us remark that Km = K. 

Corollary 5.4. Let a G (0,2]. There exists a modification X* of X^ on K such that for all j = 
1, . . . ,m 

\X*^{x) - Xl{y)\ 

lim sup , ^ , , ,^ , =0 a.s. 

no ,,yeK, \\x - y\\H^\og\\x - y\\\HAPj-^)+(^+U2+e 

\\x-y\\<S 

for any e > 0, where pj = max 1^, (3 = 1/ a if a ^ 2 and /? = if a = 2. 

l<k<j 

Proof. It follows from Theorem 15.11 and Corollary 13.31 since aj < for any j = 1 . . . d. □ 

Corollary 5.5. Let a G (0,2]. The random field X* has locally H -Holder sample paths on for 
every H G (0,i?m). 

Now let us give the directional and global Holder critical exponents of X^- 

Proposition 5.6. The random field X^ admits Hm as Holder critical exponent. 

Moreover, for any j = 1,... ,m, for any direction u G Wj, the field X^ admits Hj as directional 

regularity in the direction u. 

Proof. For Z a real SaS random variable we let 

||Z|U = (-log(E(exp(iZ))))^/". 

Then, for any x, y G M"^, 

||X*(x) - X*(y)|U = C^{iE{x - y))TE{x - y) 

where for aW. G Se 

Ca{0) = (co, [ e^<^'«> - 1 " V(C)~°~^dcl ^ and c„ = — T |cos(t)|°dt. 
V JRd J 27r Jq 

From Lebesgue's Theorem, the function Cq is continuous on the compact set Se with positive values. 

Let us denote nia = m.mCa{d) > 0. According to (fTOll in Corollary 13. 4^ for any j = 1, . . . ,m and u 

S&Se 

a direction in Wj , 

||X*(tn) - X*^{su)\\a > m^ci \t - s\"^ |log \t - s\\-^^^-'^^^ . 
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Therefore, for any s > Hj, it implies that is almost surely unbounded as |i — s| i 

so ()2ip fails almost surely on K Ci Mu. 

Moreover, Corollary I5.4l implies that {X^{tu))f-^^ satisfies (f2T]) on KnWu for any s < Hj and thus Hj 
is the directional regularity of Xa in the direction u. 

Moreover, one can find a direction u E S^~^ in which almost surely {X'^(tu))^^^ does not satisfy ()2ip 
on K n Mti for any s > H„i- Therefore, almost surely {X'^{x)) ^^-^d can not satisfy (|2T|) on K for 
any s > H^- Then, by Corollary 15.51 admits Hm as Holder critical exponent. 

□ 

Remark 5.1. In the diagonalizable case (see Example 12. 2p . the Wj are the eigenspaces associated 
with the eigenvalues of E. In particular, for 9j an eigenvector associated with the eigenvalue \j = aj, 
the critical Holder exponent in direction 6j is Hj = 1/aj. 

Proposition 15.61 compared to Theorem 5.4 of shows that the operator scaling stable field, defined 
through an harmonizable representation share the same sample paths properties as the Gaussian ones. 
Therefore it is natural to have also the same result of Theorem 5.6 [U] for the box- and the Hausdorff- 
dimensions of their graphs on a compact set. We also refer to Falconer [U] for the definitions and 
properties of box- and the Hausdorff-dimension and keep the notations of [S]. More precisely, we fix 
a compact set K CW^ and consider Q{X*){io) = {{x, X*{x){ijj))]x G K} the graph of a realization of 
the field X* over the compact K. We denote dim7^^(X*), resp dimsG^X*), the Hausdorff-dimension 
and the box-dimension of G{X^), respectively. 

Proposition 5.7. Almost surely 

dimn G{X*^) = dims g{X*) = d+l-Hm. 

Proof. The proof is very similar to those of Theorem 5.6 [6]. It also uses same kinds of arguments as 
in [3]. For sake of completeness we recall the main ideas. Corollarv 15.51 allows as usual to state the 
upper bound 

dimng{X*J < ~dI^g{X*J <d+l-Hm, a.s. 



where dimg denotes the upper box-dimension. The lower bound will also follows from Frostman 
criterion (Theorem 4.13 (a) in |9]). One has to prove that the integral Is 

Is= f ¥.\{{Xl{x)-X*^{y)f + \\x-yfy''^] dxdy, 

is finite to get that almost surely dim'^^(X*) > s. In our case, the fundamental lemma of [3] allows 
us to write this integral using the characteristic function of the SaS field X*. Actually, when one 
remarks that, using Fourier-inversion, (^^ -|- = ^ J^e"^^ fs{t)dt, where fs £ L°°(M) n 

one gets 

^s= h^lk-yir' / e 1*1 fs{t)dt]dxdy. 



KxK 



2ir' 



By a change of variables, as fs G L°°(]R), one can find C > such that 

Is < C [ \\x-y\\^-'\mx)-X:{y)\\-'dxdy 



KxK 



< 



J KxK 



where v72q, — min ^Cfj J'^d |e^^^'^^ — l| V^(^) ^ ^^^^ as introduced in the proof of Proposition 15.6 
Since J^y,j^ \\x — y\\^~^TE{x — y)~^dxdy is proved to be finite as soon as s < d -|- 1 — Hm in [6], 

dims g(X;) > dim7^C?(X*) >d+l-Hm a.s. 

and the proof is complete. 



□ 
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Harmonizable operator scaling stable random fields share many properties with Gaussian operator 
random fields. In particular, they have locally Holder sample paths and critical directional Holder 
exponent depending on the directions. In the Gaussian case (a = 2), [6] establishes such properties in 
the framework of harmonizable and moving average Gaussian operator scaling random field. However, 
for stable laws, harmonizable and moving average representations do not have the same behavior as 
we see in the next section. 



6. Moving average representation 

Let us recall the definition of moving average operator scaling stable random fields introduced 
in [6]. Let < a < 2. We consider Ma{dy) an independently scattered SaS random measure on M"^ 
with Lebesgue control measure, see [27] for details on such random measures. Let us recall that, as 
before, q = trace(£'). Let 99 : ^ [0, 00) be a continuous S-homogeneous function. We assume 
moreover that there exists (3 > 1 such that 93 is {(3, E)- admissible. According to Definition 2.7 of [6] 
it means that ip{x) 7^ for x 7^ and that for any Q < A < B there exists a constant C > such 
that, for A < \\y\\ < B, 

(22) \^{x + y)-^iy)\<CTEix)^ 

holds for any te{x) < 1. 



Definition 6.1. Since f3 > 1, the a-stahle random field 

(23) Z^{x) = [ Ux - y)!-"/- - (/^C-y)^-"/") M^{dy) , x G M". 

is well defined and called moving average operator scaling stable random field. 

Prom Theorem 3.1 and Corollary 3.2 of [6], it is stochastically continuous, has stationary increments 
and satisfies the following operator scaling property 

Vc > 0, \^Zaic^x);x e M^} ^=^^ |cZ„(rE);x G M^}, 

as the harmonizable field Xa. 

In the Gaussian case (a = 2), the variogramme of Z2 is similar to the one of the harmonizable 
field X2. Then, [6j proves that Z2 and X2 admit the same critical Holder sample paths properties. 
However, when a G (0,2), let us recall that moving average self-similar a-stable random motions 
does not have in general continuous sample paths (see [27j). The next proposition states the same 
property for Za. 

Proposition 6.1. Assume a G (0,2) and d > 2. Then, the random field Z^ is almost surely 
unbounded on every open ball. 

d 

Proof. Let us remark that 99(0) =0 by continuity and i?- homogeneity and q = ^ Oj > d > a, as 
soon as d > 2. Then, for any U open set, since f/* = C/ n Q'^ is dense in [/, for any y & U 



nU*,y)= sup <^(x-2/)i-5/"-¥,(-y)i-5/" 



+00. 



Then f* (U* ,y)'^dy = +00 and the necessary condition for sample boundedness (10.2.14) of 
Theorem 10.2.3 p.450 of [27j fails. We conclude the proof by Corollary 10.2.4 of [27j. □ 
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7. Appendix 

Proof of lemma [37^ The lower bound is straightforward. Actually, for any t > 0, is an eigenvalue 
of the matrix t'^ and therefore = \t^\ < Wt"^]]- 

Let us prove the upper bound. First, let us assume that J is a Jordan cell matrix of size I. In this 
case A = a G R and 

/I ... 0\ 
logt 1 : 
: ■•• ■•. 

V hHyr ••• log* 1 / 

Let us recall that the norm defined for a matrix A = (aij)i<ij',<(i by ||j4||oo = max X]f=i 

~ ' '~ l<i<d ■' 

subordinated norm of the infinite norm H^^Ho^ — max \ for x G M*^. By definition, we can deduce 

l<i<d 



that \\t' 



-. Then, 



\t'\\<vi\\t'\\^<vinog{t)t^Y.- 



for any t £ (0, e~^] U [e, +00). Therefore, for any t G (0, e^-*^] U [e, +00), we have 

p^ll < Vlet''\logt\^"K 

In the second case, let us assume that J is a block of the form ([7]) of size 21 associated with the 
eigenvalue A = a + z6 for 6 / 0. Then t'^ = t°-R{t)N{t) where 



/ Rbit) 



R{t) 



\ 



Rbit) 



V 



with Rhit) 



cos(Mogt) — sin(61ogt) 
sin(61ogt) cos(61ogt) 



and 



/ h 



Rb{t) j 







\ 



with Nj{t) 



V A^;-i(t) ... A^i(t) h I 



Hence, 



\t'\\<nRm\\^m- 



Since is an orthogonal matrix, ||ii(t)|| = 1. Furthermore, A^(t) is a (2/) x (2/) matrix and 



Therefore, we also obtain that 
for any t G (0, e"^] U [e, +00). 



It-^ll < V2/et"|loet|'-^ 



□ 

Proof of Proposition \3.4i Let r G (0,1). One can find rE G (0, r) such that for any ||2;|| < rE we 
have te{x) < . Let x G ©i=jQ ^fc with ||x|| < rE- Then x = te{x)^ (.e{x) and Ie{x) G ©i=j(, ^fc- 
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We first establish the lower bound of Proposition 13.41 Let us write Ie{x) = YlT=i ^fc(^) where each 
hix) G Wk- Let Lk be the coordinates of lk{x) in the basis ^/^fe-if,^;^, • • • i fj^^ of ^fc- Hence, 



by definition of P, 



P-^Ie(x) 



and X = te{x) Ie{x) = P 



/ te{xY'Li 

V TEixY'^L^, 



Since Ie{x) G 0i=,g Wt, Lfc = for /c ^ {jo, • • • 



1/2 



1/2 



< 



i^'ii Elk^(^ 



By Lemma 13.2 



1/2 



|x|| < V2e||P|| I lkTE{xf'"=\\ogTE{x)\ 

\k=jo 



2fe-l)||^^^|,2 



since te{x) < 1/e. Hence, since te{x) < 1, |logr£;(x)| > 1, > and < pj^j = max < d, 

jo<i<j 

\\x\\ < V2de||P||T£;(x)'^^o|logrs(x)|(f«--i)(n=,oll^fcf^'^' 
< v^e||P||Tij(x)"^o|logrij(x)|(P^o.i-i)||p-i/^(a;)|| 



Then, 
(24) 



< V2deME\\P\\\\P~^\\TEixp'>\logTE{x)\^^'o,j~^) 



where Me is defined by ([6]). 

Let us take the logarithm of this equation. Choosing tq < min(l,r£;) small enough, one can find 
C > such that 

(25) I logr£;(2;)| < C\ log ||x||| for ||2;|| < tq. 

Using in (j24p . we obtain that there exists C > such that for ||x|| < tq 

||x||^^o|log||x||r^^ofeo,j-i) < CTEix). 

By continuity of the map, 

X ^ ||2;||^^o|log \\x\\\~^'o{Pjo.3-^)T-^(^x)~'^ 

on {x G < ||x|| < l}, up to change C, the previous inequality holds for ||x|| < r, which gives 

the lower bound in (jlip . 

Let us now establish the upper bound in (fTTj) . We write x = J2T=i with each x^ G W^- We then 
denote by Xk the coordinates of Xk in the basis I /y-fe-i ,- , , , . . . , /y-fc r ) of Wk- Since x G W^fe, 
Xfc = for every k ^ {jo, . . . , j}. Hence, by definition of P, 

and lEix) = te{x)~^x = P 



p-'x 



( te{x)-''Xi \ 

\ TE{xy-^"^Xm J 



1/2 



Then, \\Ie{x)\\ < \\P\\ ( ^ \\TE{xr^'Xk 
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The Lemma 13.21 yields 



1/2 



Ie{x)\\ < V2e\\P\\ lkrE{x)-''''\logTE{x)\^^''-'^\\Xk 



U=30 



since r£;(x) ^ > e. Hence, since T£;(x) ^ > e > 1, < and Ik < Pjo,j^ 



0<mE < V2de||P||||p-i||rs(x)-"^|logr£;(x)r^o..-i^/^^.J|Xfc|| 



< V2de\\P\\\\p-^\\TE{x)-''^logTE{x)f''>-'~^\\P~^x\\. 
Then, using ()25p and < ||P~^||||x||, there exists a constant C > such that 

te{x) < C||x||-^^|log||x|||^^'(*'^o..-i) 
for ||x|| < ro- By continuity of the map 

^ ^ te{x) 

||x||^-'|log||x|||^-'(^^o-^~^) 

on |x G M'^/O < ||x|| < l}, up to change C, the previous inequahty holds for ||3;|| < r, which gives 
the upper bound in (fTTIl and concludes the proof. 

□ 



Proof of Lemma \5.S\ . It is sufficient to consider 

I{h) = E^m(^„)^^/"min [Me 

By definition, 



m 



m(e)^"^/"V(0"^^^'^" min ( Me 



,2 di. 



Using the formula of integration in polar coordinates with respect to E'*, see Proposition 13. H 



m 



S^t Jo 



+ 00 



m r 



^j(r^ 9) min M^; 



{hr)^*e ,2) r"~^draEt{d9). 



Since ^p is i?*-homogeneous, 



S^t Jo 



+00 



m\ r 



1-2/a 



-2-2q/a 



min ( Me 



{hrf'e ,2) r~'^+'i^^~^'''''drcjEt{d6) 



1-2/" 



/ / 7/;(^)-2-2'?/"niinfME (hrfo , 2 j r-3|log(r)P+'')(2/"-i)rf^^^,(rf0). 
J Sj^t J 



By the change of variable p = hr, I{h) is equal to 
n+00 

-2-2"/° min ( Me 



:i-2/"/i2 



5^* ^0 



p 



,2 p- 



{l+^;y){2/a-l) 



draEt{d6). 



For any 7 G (0, 1), there exists such that for every p > and every /i < 1 — 7, 

log = |log(p) - log(/i)| < A^||log(/>)| + l||log(/i)|. 



Since 2/a > 1, 
with 



/(/i)<<"-^cl-2/°/i2|log(/i)|(l+'')(2/a-l)(^^^^^) 

/i = 4 / ^(0)-2-2''/"a^.(d0) /^"p-^^llloglp)! + l|(i+^)(2/-i) 
JS^t Jl 
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and 



-2-2g/a 



aEt{de) 











f 




Jo 





p-3||log(p)| + 1 



(l+^)(2/a-l) 



dp, 



where Me and M^t are defined by ([6]). Since ^p is continuous with positive value on the compact 
set Se*, 

Hence Ii < +oo. 

It follows from Proposition 13.41 that for any 5' G (0, 1), there exists a constant c'g > such that 



II < C/i|lo, 



lZ-l 



for all p < 5' . Hence, since ai > 1, 















Jo 





p-3||log (p)| + l 



{l+r?)(2/a-l) 



dp < +00 



and I2 < +00, which concludes the proof. 



□ 
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